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The Anderson model on the Bethe lattice is historically among the first for which an 
energy regime of extended states and a separate regime of localized states could be 
established. In this paper, we review recently discovered surprises in the phase diagram. 
Among them is that even at weak disorder, the regime of diffusive transport extends well 
beyond energies of the unperturbed model into the Lifshitz tails. As will be explained, 
the mechanism for the appearance of extended states in this non-perturbative regime 
are disorder-induced resonances. We also present remaining questions concerning the 
structure of the eigenfunctions and the associated spectral statistics problem on the 
Bethe lattice. 
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1. A primer to the Anderson model 

Motivated by the quest for a theory of quantum transport in disordered media, in 
1958 AndersorP^ came up with a model for a quantum particle in a random energy 
landscape. In its simplest form, the configuration space of the particle is chosen to 
be Z d and the energy landscape is modeled by independent identically distributed 
(i.i.d.) random variables {w(x) | x € Z d }. The particles' kinetic and potential energy 
is described by the operator H\(uj) = T + A V(lu) acting in the Hilbert space £ 2 (Z d ), 
where 

• the first part is the Laplacian without its diagonal terms, i.e. (Tip)(x) := 

-T,\ v - ai \=ii>(y)> 

• the second part stands for the multiplication operator corresponding to the 
random variables, i.e. (V(uj)?p)(x) := uj{x) ij)(x), and 

• the parameter A > was introduced in order to tune the strength of the 
disorder. 

Among the interesting features of this operator is a conjectured dimension- 
dependent, energetically sharp transition from a regime of localized states to one 
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Fig. 1. Qualitative sketch of the expected phase diagram of the Anderson model for d > 2 in case 
Pa is a uniform distribution in [—1,1]. The outer straight lines mark the edges of the spectrum 
[—2d — A, 2d + A]. The inner curve indicates the conjectured mobility edge separating localized 
states from dclocalizcd ones for d > 2. 



of diffusive transport, cf. Figure [T] below. In order to describe the details of this 
transition, let us first briefly recall some mathematical framework. 

Since H\{oj) is almost surely (a.s.) self-adjoint, spectral theory is at hand for the 
analysis of the quantum time evolution e~ itH *( u ) , i n this context, it is important to 
note that the ergodicity of the random process implies that the spectrum of H\(ui), 
as well as its components in the Lebesgue decomposition are a.s. non-random, i.e. 
there is T,f C K such that for almost every (a.e.) u one has er# (7Ja(u;)) = with 
# standing for the spectrum, respectively the pure point (pp), singular continuous 
(sc) and absolutely continuous (ac) component. While it is rather straightforward 
to prove that the spectrum equals T,\ = [—2d, 2d] + A supp Pq, with Pq denoting the 
law of a single variable, determining the spectral components is only trivial in the 
limiting cases: 

• the Laplacian has purely ac spectrum in a(T) = [—2d, 2d] and its general- 
ized eigenfunctions are extended plane waves. 

• the random multiplication operator has a.s. only (dense) pp spectrum, 
a(V(u))) — {oj{x) | x e Z d } = suppP an d its eigenfunctions {5 X | x G Z d } 
are localized at the lattice sites. 

Within the framework of spectral theory, one may now speak of spectral localiza- 
tion respectively derealization within an energy interval / C R if the operator 
H\(uj) has a.s. only pp spectrum respectively only ac (or sc) spectrum there. While 
the RAGE theorem supports some identification of the spectral with dynamical 
properties, stronger dynamical notions of localization respectively derealization ex- 
ist. For example, one speaks of strong exponential dynamical localization within 
I C M if initially localized wave packets with energies in / remain localized: 



\x\=R 



mp\(6 x ,e- ita »Pi(Hx)6o) 



<Ce- R /^ (1) 
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where £/ > is the energy dependent localization length. Transport on the other 
hand is often captured by a transport index a € (0, 1] which describes the asymp- 
totic broadening of initially localized wave packets with energies in I C K in the 
sense that e.g. the second moment 

(\x(t)\ 2 ) := ^ Ix^Efl^^e-^Pz^)^)! 2 ] (2) 

grows asymtotically as Cf 2a for t — >• oo. Ballistic transport as present in the case 
A = corresponds to a = 1. In case A > 0, one conjectures that the transport, if 
present at all, is always diffusive, i.e. a = 1/2. 

Many other expressions of transport such as a non- vanishing conductivity tensor 
exist and are investigated. We refer the interested reader to the textbooks' 19 ! 40 ! 44 ! 
and further articles^!!] on the subject. 

Returning to the conjectured spectral and dynamical properties of the operator 
H\{uS), the situation is summarized as follows: 

• for d = 1, 2 one expects spectral (and dynamical) localization throughout 
the entire spectrum for any A > 0. 

• for d > 2 the situation is best summarized by the phase diagram in Figure[T] 
Its essential feature is a sharp separation of a regime of spectral and strong 
exponential dynamical localization from that of diffusive transport by a 
so-called mobility edge. 

Localization in low dimension (d = 1,2) is clearly a non-perturbative quantum 
phenomenon. While the behavior for d = 1 has been established by means of transfer 
matrix methods (starting with the works by Goldsheid, Molchanov and PastmP^), 
establishing complete spectral localization for small A > in dimension d = 2 
remains one of the open problems in the field. 

In higher dimensions (d > 2) localized states occur at the edges of the spectrum, 
where the states arise from large deviations of the potential, and in the region of 
large disorder in which the features of the multiplication operator dominate. These 
regimes were made amenable to a mathematical proof through the celebrated multi- 
scale analysis by Frohlich and Spencei^who first showed the absence of ac spectrum 
in these regimes. Based on this result, the proof of pp spectrum was later completed 
by the analysis of Simon and Wolff!^ A proof of strong exponential localization at 
band edges and large disorder was presented by Aizenman in Ref. 5, It uses the 
fractional moment method developed by Aizenman and Molchanov!^! 

Generally one expects that some ac spectrum remains at small disorder in case 
d > 2. As already mentioned, these extended states should go along with diffusive 
transport. Establishing any of these properties remains a huge challenge in the field. 
Steps towards an understanding have been made by Erdos, Salmhofer and Yau, who 
proved diffusive behavior in a scaling limit P^l 
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Let us close this short review by briefly mentioning two works related to the 
quest of an understanding of derealization of quantum particles in random media. 
Disertori, Spencer and Zirnbauer have successfully established the existence of a 
diffusive phase in the supersymmetric hyperbolic non-linear sigma-modelP^ Very 
recently Erdos, Knowles, Yau and Yin established derealization and a diffusion 
profile for the eigenvectors of certain random band matrices.^ 



2. The phase digram on the Bethe lattice 
2.1. Some history and the question 

In view of the difficulties of establishing the existence of a regime of derealization 
for any high dimensional lattice Z d , it is only natural to look for simpler models. 
Historically among the first was the study of the operator H\(ui) = T + AV(w) 
in the Hilbert space £ 2 (M) over the so-called Bethe lattice. The latter is a regular 
tree graph, i.e. a graph without loops in which every site has the same number 
K + 1 ( > 3) of neighbors. The almost-sure spectrum of the operator H\(ui) on this 
graph is easily calculated: 

S A = [-2y/K,2y/K] + AsuppPo- (3) 

It was Abou-Chacra, Anderson and Thouless who first looked at this model 



and proposed the so-called self-consistency equations (cf. (20) below) as an ap- 
proximation to behavior of the self-energy in high dimensions lil^l Analyzing these 
self-consistency equations, they presented convincing arguments for 

• the presistence of a regime of delocalized states for small disorder within 
the ac spectrum of the Laplacian, er(T) = [—2>/K, 2\J~K\. 

• the occurrence of localized states for large disorder and outside the interval 
[—(K + 1), (K + 1)] for small disorder. In view of the latter is a mean- 
ingfull statement only if suppPo = i.e. for unbounded random variables. 

Both of these statements were later proven: the localization result in the strong sense 
of ([I]) by AizenmarP and the persistence of ac spectrum alongside ballistic transport 
by KleinlHEil (Different methods establishing the persistence of ac spectrum can be 
found in Refs. |7p5| ) 



Both of these results left open the question of the location of the mobility gap 
for small disorder. In particular in case supp Pq = R, the nature of the states in the 
energy regime 2\J~K < \E\ < K + 1 remained undetermined for small A. This ques- 
tion was already taken up by Miller and Derridap^ but their (mostly numerical) 
results remained inconclusive. The situation is summarized in Figure [2] 

Triggered by interests in many-particle localization for which the tree is proposed 
as emulating features of the configuration space of many particles,^ the question of 
the location of the mobility edge has recently received renewed interest in the physics 
community. For example, the phase diagram in case Pq is a uniform distribution on 
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Fig. 2. A sketch of the previously known parts of the phase diagram. The outer region is of proven 
localization in the sense of only pp spectrum and strong exponential dynamical localization. The 
smaller hatched region is of proven derealization in the sense of purely ac spectrum and ballistic 
(= diffusive, cf. Subsection \2.3\ transport. The new results address the question of the nature of 
the states in the middle. 



[—1, 1] was numerically investigated in Ref. 14 which presented a diagram resembling 
the features of the one on the regular lattice (cf. Figure [3]) . Before we come to a 
critique of the results of Ref. [14] and the solution to the riddle of the location of the 
mobility edge for small disorder, let us first make a detour into other related results 
in the literature which present some red herrings. 

In case of unbounded random variables, supp P = M, the states in the regime of 
small disorder A and energies \E\ > 2\/~K stem for large and hence rare fluctuations 
of the random variables. As a consequence, the density of states D X (E), i.e. the 
Lebesgue density of the integrated density of statesQ 

N X (E) := E[(8 ,P ( _ OOiE) (H x )6 )]= f D x (E')dE' , (4) 

J — oo 

is expected to be severely suppressed. This phenomenon is usually referred to as 
Lifshitz tailing. In fact, it is not hard to show that in case of Gaussian random 
variables, N X {E) < exp (-C(E)/X 2 ) for any E < -2sfK and all A > OP Similarly, 
in case of bounded non-negative random variables ui(x) > 0, it is conjectured thalP^ 

Um loglog|logjVA(-2^ + £)| = _ I 

eio log E 2 ' ^ ' 

Since the rare fluctuation which cause the suppression of the density of states are 
typically spatially separated, it is quite natural to suppose that the states in this 



a Note that the existence of a Lebesgue density is guaranteed by a Wegner estimate in case Pq 
has a bounded density|33ln fact, if Pq is sufficiently close to a Cauchy distribution, the density of 
states in known to be analytical 
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energy regime are localized in isolated wells. As we will point out next, this is not 
the case and there is a regime of ac spectrum well within the regimes of Lifshitz 
tails. 



2.2. The location of the mobility edge 

The phase diagram of the Anderson model on the Bethe lattice turns out to be 
described in terms of the free-energy function 

logE[\G x (0,x;z)\ s ] 
ipx[s;z) := hm r-j , (6) 

|x|— yoo \X\ 

which captures the decay properties of the Green function G\(x, y; z) :— (S x , (H\ — 
z)~ 1 8 y ). Based on factorization properties of the Green function, it is proven in 
Ref.[9]that this free energy function is well defined, monotone decreasing and convex 
in s > for all z G C + := {w G C\lmw > 0}. The following properties and 
quantities are essential for the description of the spectral phase diagram: 

(i) For s G (0, 1) the limits \x\ — > oo and Imz 4- exist and commute. The boundary 
limit 

<P\(1;E) =lim]imcp x (s;E + ir)) (7) 

sfl rjj.0 

exists for Lebesgue almost all E G K. 

(ii) The derivative at s = can be identified with a Lyapunov exponent 

L x [z) := -lim^M = _Epog|r A (0;*)|] , (8) 

sio s 

where the last equality results from the factorization of the Green function, and 
T\(0;z) :— (6o,(Hf — z)~ 1 5o) denotes the Green function of the operator H\ 
restricted to the Hilbert space over a regular rooted tree T with root denoted by 
0. 

For the formulation of the result, which was proven in Ref. [9j let us briefly state 
the assumption on the probability distribution Pq which will be assumed for the 
validity of all of the theorems in this review: 

Assumptions on Pq: The single-site probability measure has a bounded density, 
i.e. Po(dv) = g(v)dv for some g G L°°(M), which satisfies: 

(i) it is bounded relative to its minimal function, i.e., for Lebesgue-a.e. »€l 
and some c < oo: 

g(v) < c inf (2i/) _1 / l\ x - v \< v g(x) dx . (9) 

"6(0,1] J 

(ii) the moment condition with some r > 0: 

\v\ r g{v)dv < oo. (10) 
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The first theorem determines the location of the mobility edge in terms of 
<Px(l;E). For purely technical reasons this theorem was proven in Ref. 9 only in 
the unbounded case. The second part of the theorem holds more generally and was 
already established in Ref. [5j 

Theorem 2.1 (Characterization of the mobility edge). If infi^i<jt q(v) > 
for all k < oo, the free-energy function determines the mobility edge in the sense 
that 

• for a.e. E £ Y,\ ifip.\(l;E) > — logK then a.s. 

DI c (E,uj) := -lmG x (0,0;E +iO,u) > 0. (11) 

• if Va(1; E) < — logK for a.e. E € I C then I is a regime of strong 
exponential dynamical localization in the sense of (JlJ. 



The left-hand side of (11) is the density of ac states at energy E. The event 
D^ C (E, lj) > happens with either full probability or probability one. Moreover, by 
translation invariance, if D^ c (E,uj) > a.s., then the ac density associated to any 
other sites is stricly positive, Im G\(x, x; E + zO, u) > a.s. 

The above theorem identifies the set {(E, A) G R 2 |^a(1;£0 = — log if} as the 
location of the mobility edge for the Anderson model on the Bethe lattice and hence 
paints a fairly complete picture. Unfortunately, due to a lack of regularity results 
the question whether this set constitutes a line in the phase diagram remains open. 
In fact, since the function tp\(l;E) is hard to analyze, we do not employ the above 
criterion directly to answer the question about the location of the mobility edge 
for small A. However, in this regime it is reasonable to expect that the estimate 
(px(l;E) > —L\(E), which derives from the convexity of the free energy function, 
is rather tight. The following theorem (which for unbounded random variables may 
also be viewed as a corollary to Theorem |2.1| but under stronger hypothesis is 
proven m Ref. 19} summarizes the resulting criterion for extended states. 

Theorem 2.2 (Lyapunov exponent criterion). For a.e. EE Y,\ the condition 
L X (E) < log K implies D a x c (E) > a.s. 

Based on this theorem, the determination of the mobility edge for small A now 
proceeds along the following lines: 



(i) In case A = the Lyapunov exponent can be explicitly calculated (using ( 20 1 
below). One obtains: 

(=logVK ii\E\<2VK, 
£ (log VK, log K^j $2VK<\E\<K+1, (12) 
> log K if \E\ > K + l. 
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As an aside we note that since Go(0,x;z) = Go(0, 0; z) e> x > L °( z \ the explicit 
expression helps identifying E = ±(A + 1) as the threshold for the summability 
of the unperturbed Green function. In other words, the interval [— (K+l), K+l] 
is the spectrum of T on the Banach space ^ 1 (B). 
(ii) Using the fact that the Lyapunov exponent is the negative real part of a Herglotz 
function, it is not hard to convince oneself that it is weakly continuous in A in 
the sense that for any interval Id: 

lim / L x {E)dE= [ L (E) dE . (13) 
-Uo J i J i 



In view ( 12 ) this implies^ the existence of some ac spectrum within the 



full interval [—(K + 1), K + 1) for small A and supp g = K. 

In case of Cauchy random variables, i.e. g(v) — 7r _1 (i> 2 + one can even 

make a stronger statement since the Lyapunov exponent can be explicitly com- 
puted: L\(E) = Lq(E + iX). The Lyapunov exponent criterion hence identifies 
a full interval within (—{K + 1),K + 1) as the regime of ac states for small 
A. Since in this case the ac density is strictly positive for a.e. energy in a full 
interval, by a Simon- Wollf type argument, one may even conclude that the ac 
spectrum on this interval is pureP 
(iii) In case of bounded random variables, e.g. supp q = [— 1, 1], we prove^the conti- 
nuity of the Lyapunov exponent at the spectral edge E\ := inf Y,\ = — 2sfK — A: 

limsup£ A (£0 < L {E X -X). (14) 

As a consequence, for any < A < A m , n := (V K — l) 2 /2 there is some 5\ > 
such that the Lyapunov condition is satisfied on the full interval [E\, E\ + 6\] 
and hence there is purely ac spectrum occurring at the band edges within 
the regime of conjectured Lifshitz tails |5]). In particular, this shows that the 



numerical results in Ref. 14 need to be corrected for small A. In view of (131 or 
the perturbative results of Klein p2l which cover energies within [—2>J~K, 2\fK], 
we conjecture that for sufficiently small A the entire spectrum is purely ac, 
cf. Fig. |3j 

It is of some interest to estimate the critical disorder strength A c > above 
which the entire spectrum is pp and dynamical localization prevails. Using the 
factorization property of the Green function together with bounds of the form 

E[\G x (0,0;z)\ s \u(x),xjtO] < sup / . do < J |g|l °V =: C S (X) 



\Xv - a\ s ~ (1 - s) X s 



on the conditional fractional moment of the diagonal of Green functions, conditioned 
on all random variables aside from u)(0), one obtains (<f(l; z) < ) ip(s; z) < log C S (A). 
Optimizing over s € (0, 1), i.e. supposing that A > e^Hoo and picking (1 — s) _1 = 
log(A/||e||oo), we arrive at 

A c < HellooA^elogA'+l) . (15) 
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Fig. 3. The solid line sketched the previously expected phase diagram for the Anderson model on 
the Bethe lattice (cf. Ref |14| in case Pq is a uniform distribution in [—1,1]. Our analysis suggests 
that at weak disorder there is no localization and the spectrum is purely ac. While the proof of 
that is incomplete, we do prove that for A < {\/K — l) 2 /2 there is no mobility edge beyond which 
localization sets in. 

This bound is consistent with the conjectured asymptotics A c ~ CKlogK for 
large K. 

A lower bound on A c is only available in case of Cauchy random variables (or, 
more generally, if the random variables have a Cauchy component). It is based 
on the Lyapunov exponent criterion and the fact that L\(E) = Lq(E + iX). The 
equality Lo(i\) — log if yields A c > K — 1 in this case. 



2.3. Transport properties 

In view of the existence of a regime of extended states within the regime of extremely 
low density of states, one might wonder whether there is a further dynamical tran- 
sition if one crosses from the perturbative regime already analyzed by KleirP^ into 
that of Lifshitz tails. The following theorem excludes this possibility and establishes 
ballistic transport throughout the region of positive ac density. 



Theorem 2.3 (Ballistic transport in the ac regime). Assume that (10) holds 
with r = 12. Then for any initial state tp = f(H)S Q , with f G L 2 (R) supported in 
{ E e K | D a x c (E) > a.s.} and all b > 0: 



E 



\x\<b/r, 



dt 



< C(f)b + o{n) 



(16) 



with some C(/) < oo, and o(rj) a quantity which vanishes for n — > 0. 
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Theorem |2.3| in particular implies that within the regime of ac states the time 
evolution is at least ballistic in the time averaged sense 

p oo 

liminfr/ 3 / e^ 1 (\x(t)\ 2 ) dt > . (17) 

This proves the ballistic nature of the ac states, since quite generally quantum 
evolution is at most ballistic, limsup^g rf J^°e~ vt (\x(t)\ 2 ) dt < oo, . 

One might wonder whether this result fits the general conjecture about the 
diffusive nature of the transport in random media. For this it is useful to note that 
on the Bethe lattice the diffusion is ballistic: 

(6 x ,e tA S )dt=(6 x ,(-A)- 1 6 a ) = w (18) 

where (Aip)(x) := —Tip(y) — (K + l)ip(x) stands for the diffusion generator. The 
proof of Theorem |2.3| which can be found in Ref. '12} essentially proceeds through 
the following estimate on the second moment of Green function 

sup E[|G A (0,x;C)| 2 ] (19) 



for any / C {E e R\D<l c (E) > a.s.}. In view of ([18]) this bound reflects the 
diffusive behavior. 



3. Resonant derealization 

We will not attempt to give a proof of any of the results presented in the previous 
section, but only sketch the ideas behind the geometric resonance mechanism which 
we dubbed resonant derealization and which is responsible for the non-perturbative 
appearance of ac states outside c(T). To simplify the presentation, we restrict ourself 
to the case of a regular rooted tree T with root 0, for which every site has K forward 
neighbors. In this set up, the truncated resolvents T\(u;z) := (5 U , (H^ u — z)~ 1 8 u ) 
associated with the subtree T+ which is rooted at and forward to u satisfy the 
so-called self-consistency equations: 

Tx{x;z)= (\w{x)-z- Vr A ( K ^ . (20) 



Here H x u denotes the restriction of H\ to £ 2 (T+). The equations (20) consitute 



a discrete dynamical system system for probability distribution on C + induced by 
the map u> i— > r^(0; z,lu). In order to establish the presence ac spectrum at energy 



E within the spectrum, it hence remains to investigate the fixed points of ( 20 1 for 
the boundary values z = E + iO and establish the positivity of the right-hand side 
of 

Imr A (0;z)> \Gx{0,x;z)\ 2 lmT x (x + ;z), x ~< x+ . (21) 

\x+\=R 
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Here the sum ranges over the sphere of sites at distance R from the root in T and 
x -< x + denotes the unique predecessor of x + . 

We will now present an argument, which proves the following implication: 

If lmr,\(0, E + iO, uj) = a.s. and L\(E) < log -FT, then there exists some S > 
such that the event 

Nr(uj) 1= ^ l {Gx(0,x;E+iO,uj)> e*M} > 1 ( 22 ) 
\x\=R 

occurs with a positive probability, independently how large R > is. 

First and foremost, this shows that the Green function cannot be square summablc 
(and hence by a Simon- Wolff argument!^ there is no pp spectrum). A more detailed 



analysis based on (21) establishes the existence of ac spectrum, in contradiction 



with the assumption in the above implication. 

The idea to construct the above enormous fluctuations is based on the following 
properties of Green functions on trees: 

(i) Factorization of the Green function: G\(Q,x;z) — Ga(0,x_;z) G\(x,x; z), 
where G\ pertains to the operator restricted to the tree which is cut beyond 
x. 

(ii) Local dependence on u)(x): G\(x,x;z) = (Xoj(x) — <J\{x\ z)) , where 
(T\(x; z) is the so-called self-energy at x which is independent of uj(x). 

Since the typical decay of the Green function is given by the Lyapunov expo- 
nent, the event K x := {G x (0,x^;E + iO) > C e ~ L ^ E ^ x \} occurs asymptotically 
with probability one. In order to compensate for this typical decay, we look at the 
extreme event E x :— {G\(x : x; E + iO) > e' iA ^' +5 )l :r l} which using its local depen- 
dence on uj(x) may be proven to have a conditional probability which is of order 
e -(L\(E)+8)\x\ _ ^ g a conS equence, 

E [N R ] > n e *) ~ C K R e-( L *W+V R > 1 

|x|=fl 

provided L\(E) < \ogK and S > is sufficiently small. 

In order to establish the asymptotic occurrence of the event ( 22 ) , we use the 



second moment method which is based on the Paley, Zygmung inequality 

r ^>-»>- E M- <23) 

An upper bound on the second moment is essentially based on the following con- 
sequence of a rank-2 calculation of the conditional expectation, conditioned on all 
random variables aside from x and y: 

F(e x n£ y \u(0,Z<?{x,y}) < Cr- 1 (r- 1 +min{l ) |e(x_ J y_;^|}) . (24) 
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Here r = e^- Lx ^ E ' +s ^' x ' and G{x-, y_; E) denotes the Green function on the subtree 
of T which is chopped beyond x and y. (For more details, see Ref. |9|) 

4. Some remaining questions 

The above geometric resonance mechanics should also be applicable to homogeneous 
hyperbolic graphs other than trees. In this context, let us mention that the question 
of persistence of ac spectra in the perturbative regime has been studied for so-called 
trees of finite cone-typep^ certain decorations of the regular tre^MES] ag we ^ ag f or 
the so-called Bethe strip 1^123 It would be desirable to extend the above analysis to 
such hyperbolic graphs with loops. 

Aside from this obvious question about extcndability, there is another aspect 
to the Anderson model on the Bethe lattice which is worth clarifying and which 
starts to attract attention in the physics community) 17 ! 20 ! This concerns the so- 
called spectral statistics and the associated question of the participation ratio of 
the eigenfunctions. 

In order to put these questions into come context, let us switch back to the 
Anderson model on 1 d . Carving out a finite cube A C 1 d , one may study the random 
eigenvalues {E n (A,uj) | n e N} of the restriction of H\(uj) to £ 2 (A). The existence 
of the density of states (which is guaranteed by a Wegner estimate^) proves that 
the average density of these eigenvalues within the spectrum is proportional to the 
volume |A|. It is therefore reasonable to fix an energy E E Y,\ and look at the 
rescaled process of eigenvalues in a window of size |A| _1 , i.e. the random measure 

VA : =z2 S \A\(E n (A)-E)- (25) 

n 

The spectral statistics conjecture states that /z^ converges (weakly) as A f Z d to 

• a Poisson process if E £ E^ p . 

• a GOE process if E € X^ c . 

Here GOE process refers to the (rescaled) process of bulk eigenvalues of a random 
matrix in the Gaussian Orthogonal Ensemble. Substituting the assumption E € £ pp 
by a fractional-moment localization estimate, the first part of this conjecture was 
proven by MinamPSl for arbitrary dimensions and, for d = 1, earlier by MolchanovP^l 

The question which suggests itself is that of the spectral statistics on the Bethe 
lattice. Let me first describe how not to approach the spectral-statistics conjecture. 

If one considers a finite regular rooted tree of L generations, i.e. that subtree 
Tl C T in which all sites have distance at most L E N from the root, one may 
look at the eigenvalues {E n (L, oS) \ n € N} of the restriction of H\(oj) to £ 2 (Tl) and 
consider the random measure 

Mf : = / ] $ \r t \(E n (L)-E) ■ ( 26 ) 
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It was shown in Ref. [8] that (under a regularity assumption on the Lyapunov ex- 
ponent, which e.g. is satisfied for Cauchy random variables) the limiting process is 
always Poisson, irrespectively whether E was taken from the regime of localized or 
delocalized states. 

Theorem 4.1 (Possion statistics throughout the spectrum). Suppose that 
the Lyapunov exponent L\(E) is an equicontinuous function of E G I. Then the 
reseated process of eigenvalues fif converges weakly as L oo to a Poisson process 
for a.e. E e /. 

At first look, this seems surprising. However, as is explained in detail in Ref. [8j 
it does not contradict the spectral-statistics conjecture described above if carefully 
interpreted. The relevant limit of finite trees is not the infinite homogenous tree 
graph T, but rather a single-ended canopy graph. On this tree graph, the random 
operator H\{to) is proven to have only pp spectrum at any A > 0. 

In view of this negative result, which is caused by the presence of a large surface 
in the truncated tree, one needs to look at other finite graphs which are locally tree 
like. One example is the ensemble of random (K + l)-regular graphs. It consists 
of the uniform probability measure on graphs on N vertices where each vertex has 
(if +1) neighbors.^ As N — > oo the girth of the graph, i.e. the minimal loop length, 
diverges in probability. In fact, numerical simulations suggest that for large K the 
eigenvalue spacing distribution of the operator T approaches that of GOEpSl 

The ensemble of random regular graphs is hence well suited to investigate the 
spectral statistics conjecture of the operator H\(cu) on the Bethe lattice. The gap 
statistics of the process of eigenvalues as well as the localization properties of the 
random eigenfunctions are currently under investigation in the physics community. 
There is numerical evidence^ that the delocalized phase {{E, A) | E £ S^ c } decom- 
poses into two distinct regimes: 

• One regime (E, A) (the so-called non-ergodic phase) in which the eigenvalue 
statistics is Poisson and the eigenfunctions occupy only a fraction of the graph. 

• another regime (E, A) (the so-called ergodic phase) where the eigenvalue statis- 
tics is GOE and the eigenfunctions occupy the whole graph. 

A proof of even part of this exciting proposal remains a challenge. 
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